Asymptotic analysis of the Z-scan trace to determine 
multiphoton absorption coefficients 
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We present a theoretical analysis of open aperture Z-scan signatures for materials exhibiting an 
absorption admixture of different multiphoton processes. Based on a polynomial expansion, we 
suggest some procedures to separate the contribution to the trace of the recorded transmittance 
made by absorptions of distinct order. The signal processing from experimental data allows a 
sequential estimation of assorted nonlinear absorption coefficients accurately. 



I. INTRODUCTION 

Multiphoton absorption of intense wavefields in matter 
is a phenomenon of interest in a great variety of applica- 
tions. In multiphoton microscopy, for example, enhanced 
optical sectioning is achieved leading to superresolved 
three-dimensional imaging without the use of a confocal 
aperture. l2[ In this case, higher order nonlinear ab- 
sorption (NLA) processes are mostly proposed to further 
improve spatial resolution and light penetration. 0, 0] In 
optical limiting, on other hand, highly nonlinear absorb- 
ing materials are also expected to exhibit large transmit- 
tance at low intensities, resistance to laser-induced dam- 
age, stability over time, and ultrafast responses, [sl, 0, 0] 

In pursue of characterizing the nonlinear response of a 
material, open aperture Z-scan provides the standards to 
determine the multiphoton absorption coefficients, [s*] In 
this technique, the transmittance of a thin sample is mea- 
sured in the far field when it is scanned through the focal 
volume of a focused pulse beam. Profiting from the char- 
acteristic theoretical curve drawn at a given multiphoton 
absorption process, the corresponding NLA coefficient is 
subsequently determined by fitting the absorptive Z-scan 
trace. When three or more photons are involved in the 
NLA process, however, the open aperture scan by itself 
IS not very sensitive. 

This small sensitivity is of special significance when 
more than one nonlinear process is active. ^2] If a pulse 
beam induces a NLA effect involving a large number of 
photons, a lower-order absorption process may come out 
simultaneously within some narrow spectral windows. 
For instance, ultrashort pulse propagation at 1320 nm 
in the polydiacetylene-PTS exhibits simultaneous 2PA 
(two-photon absorption) and SPA. [13] Also, an exact de- 
generacy between three- and four-photon absorption is 
observed at the peak wavelength of 1890 nm. In chalco- 
genide glasses, it is found an important contribution of 
2PA in measuring the three-photon NLA coefficient at 
1.06 jiia in the picosecond range. [14 Il5] Even four- and 
five-photon absorption is observed in Te02 -based sam- 
ples at a wavelength of 1550 nm with a small linear ab- 
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sorption coefficient, 

In the examples given above, there is no significant dif- 
ference between the experimental Z-scan and the fits ob- 
tained using n- and (n + l)-photon absorption coefficients 
as the fitting parameter. To distinguish between these 
two absorptions, the Z-scan experiment is commonly ac- 
complished at various irradiances. [il 0, Ei, [13, El For 

sufficiently low intensities, the (n-h l)-photon absorption 
process cannot participate in the Z-scan signature; on 
the contrary this effect clearly dominates over the absorp- 
tance peak of the curve at high irradiances. Nevertheless, 
one single Z-scan trace measured at an intermediate ir- 
radiance might provide the required NLA coefficients si- 
multaneously. Pointing upon this direction, analytical or 
semi- analytical approximations are proposed to sim- 
plify this labour. 

In this paper we address the problem of identifying 
and separating the contributions of degenerate multipho- 
ton absorption processes from a single Z-scan signature. 
Our analysis exploits the fact that Z-scan traces yield 
accurate estimates of NLA coefficients if the curve fitting 
is performed with experimental data taken with sufficient 
precision. The paper is organized as follows. In Sec.HIlwe 
review the basic grounds on multiphoton absorption oc- 
curring in a thin slab made of a nonlinear material. Ana- 
lytical expressions of open aperture Z-scan signatures for 
Gaussian pulse beams are deduced in Sec. lIIIl In Sec. lIVI 
we interpret the Z-scan traces conveniently and we give 
some procedures to achieve accurate estimates for the 
order of the nonlinearity and its corresponding absorp- 
tion coefficient of the different NLA processes. Finally, 
in Sec. [Vl the main conclusions are outlined. 



II. NONLINEAR ABSORPTION BY A THIN 
MEDIUM 

Let us consider a pulsed beam propagating in air along 
the z axis and impinging normally onto a thin bulk 
medium to prove its response upon NLA. In this study 
we assume that if the beam intensity is low enough to ne- 
glect NLA in the sample, the transverse profile is not dis- 
turbed significantly due to diffraction. Exceptionally, lin- 
ear absorption in the thin sample is regarded. Addition- 
ally, we consider ultrashort pulse excitations so that slow 
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(cumulative) nonlinearities such as thermal nonlinearity 
are also ignored. [20, 21] Thus we examine m-photon ab- 
sorption and /-photon absorption simultaneously, being 
2 < m, L Starting from the slowly varying envelope ap- 
proximation we may simplify the evaluation of the field 
observed at the exit plane by means of the following dif- 
ferential equation [22] 



(1) 



and the contour condition / = 1 at Cm = 0. We point 
out that a linear absorption process may be included in 
Eq. ([1]) if / = 1 (or m = 1) and, alternatively, if we insert 
the term —rfrnil on the right-hand side of the equation, 
cases that are considered below. Here, / = I/Iq stands 
for the ratio of the output intensity / and the input in- 
tensity /o in the sample, the reduced axial coordinate 
C,m = L/Lm is the ratio of the sample depth L and the 
characteristic length for the m-photon absorption pro- 
cess. 



1 rl — m 
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being am the m-photon nonlinear absorption coeffi- 
cient. For instance, in polydiacetylenes typically [13.] 
a2 = 6 cm/GW and as = 2.7 cm^GW^ at 1320 nm 
so that L2 — 830 /im and L3 = 930 jam for a peak inten- 
sity Io = 2 GW/cm^ Finally, 
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represents the relative strength of both NLA processes. 
In our numerical example 7723 = 0.9 suggesting a balanced 
charge of both nonlinear phenomena. 

The permutation m ^ I m. ^ provides an equiva- 
lent form of the differential equation; if one combination 
shows rjmi > 1 the counterpart gives rfim = < 1. We 
may restrict our problem considering mPA a dominant 
process over /PA, whence leaving the unity as the maxi- 
mum value of the coupling parameter i^mi- According to 
Eq. (|3]), rjmi < 1 at sufficiently low intensities if m < / . 

Validity of Eq. ([1]) is restricted to nonlinear interac- 
tions with matter producing small and modest alterations 
of the beam profile. Otherwise, multiphoton absorption 
may lead to self-focusing within the sample. In mathe- 
matical terms we impose Cml^Cm^l ^ ^- Using Eqs. ([2]) 
and (|3]) we infer a sufficient condition L <C (or Cn <^ 1) 
for n = m^l. 

In general, one cannot find an analytical solution of 
Eq. ([1]). For a single NLA process being r]mi = 0, a 
solution may be provided in a closed form. In this case, 
Eq. ([1]) is reduced to d(^^I = — and the exit intensity 
in the medium yields 



(4) 



-Vl + (m-l)Cm* 



In some circumstances, e.g. (m — l)Cm > 1, Eq. (|4]) 
provides a significant beam depletion due to strong m- 
photon absorption, a result that is out of the validity in 



our approach. Therefore we focus our attention on the 
regime Cm 1- In this case we may expand Eq. (|4]) into 
Taylor series around Cm = giving 
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where 5 is a positive integer. In the limit /S ^ 00, /^^^ 
provides an exact solution of the transmitted field inten- 
sity. We point out that Cm = stands for sufficiently thin 
samples but also for sufficiently low input intensities such 
that NLA of m photons is negligible, in accordance to the 
definition Cm = (^mLI^~^. Specifically /*^^^ = 1— Cm, and 
the parabolic approximation yields 
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Cm + 'fC iVml = 0). (6) 



Interestingly, we may provide a solution of the more 
complex equation ([T]) also in terms of polynomials of Cm- 
We point out that a different approach using the method 
of Adomian may provide an approximate solution of ([T]) 
given also in terms of a series expansion. [23, 24] Develop- 
ing / into a Taylor series about the origin, the first-order 
solution yields 



/W =/o + (5c„/)oC™ = l-(l 



VmlXr, 



(7) 



where subindex stands for the specific value at Cm = 0. 
Eqs. ([5]) for 5 = 1 [also shown in (|6])] and (O are for- 
mally equivalent within this linear approach; however, 
for simultaneous NLA processes we use (1 + r]mi)Cm in- 
stead of Cm- In this context we interpret that, in a low- 
order approximation, our solution for competing NLAs 
might correspond to the solution for a single one with 
Cm L/Lmi and characteristic absorption depth 
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The average length Lmi approaches Lm in the low- 
intensity regime where Lm <C Li^ but shifts toward Li 
as long as the input intensity increases. In general, since 
two NLA processes are present, Lmi is lower than either 
Lm or Li individually. The quadratic (second-order) ap- 
proximation obtained from the Taylor series yields 



j{2) ^ f(i) , {'^ ^ r]mi){m ^ Irjmi) .2 
^ 2 



(9) 



We point out that the last term of the polynomial expan- 
sion is evaluated using the chain rule 



i = dj{dcJ)dcJ 



(10) 



and Eq. ([T]). Once again, Eq. (|9]) would be associated 
with beam depletion under a single NLA process [see (|6])] 
in using the axial coordinate (1 + r]mi)Cm^ as previously 
discussed within the linear approximation, and also sub- 
stituting the number m of photons by the average 



m = 



m + lr]mi 
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(11) 
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FIG. 1: Evolution of / when mPA and /PA are present si- 
multaneously. Depth Cm is in the horizontal axis and the 
parameter rirni = 1- Numerical evaluation of {[T]) is plotted 
in solid lines. Also, the linear approximation d?]) is drawn in 
dotted lines, and 7^^-* of Eq. ([9]) is represented in dashed lines. 



This parameter spans from m to / as rjmi moves from zero 
to infinity {r]im ^0); at the intermediate value r]rni = 1, 
fh stands for the arithmetic mean of m and /. 

Solution of Eq. ([1]) is computed numerically and rep- 
resented graphically in Fig. [T] for different values of m 
and I. A unit value of rjmi is selected in order to ana- 
lyze NLAs with balanced strengths within the sample. 
The linear approach given in ([7j) overestimates field at- 
tenuation due to absorption, whereas the quadratic ap- 
proximation (O underestimates it. The validity of these 
approximations is restricted to low depths, Cm 1, and 
deviations upon the exact solution become severe if Cm 
approaches unity. Though not represented in the figure, 
convergence with higher-order solutions results extremely 
slow showing strong oscillations upon the order S of the 
approach if Cm ^ 1- This suggests that a Fourier ex- 
pansion in terms or harmonic functions could be more 
convenient than a Taylor series of polynomials. Never- 
theless we follow a different approach next. 

Let us exploit the fact that beam depletion under si- 
multaneous mPA and /PA behaves similarly to single 
mPA, for which the exact solution (|4]) might be provided. 
Hence, using the axial coordinate (1 + r]mi)Cm instead of 
Cm in Eq- Q 5 and substituting the number m of photons 
absorbed in the process by the average fh given in (pTj), 
we finally have 



, (12) 

— {/l + (m-l)(l+7?mOCm 

In the limit r]mi = 0, obviously, Eq. (p!2|) simplifies to 
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FIG. 2: Behavior of / within depths Cm < 1, assuming mPA 

and /PA occur at T]rnl — 

1. Numerical evaluation of ([1]) is 
plotted in solid lines. Also, the analytical approach given in 
(|12p is drawn in dotted lines. 



Eq. dH) giving the exact solution. 

In Fig. [2] we inspect graphically the accuracy of (p!2|) for 
the intermediate value rjmi = 1- Within the interval of 
interest, < Cm ^ 1, deviations of the computed inten- 
sities and those available from Eq. (p!2|) have the upper 
bound 9.7 10"^ found at C2 = 1 for mixed 2PA and SPA 
processes, representing a relative error of less than 1%. 
For simultaneous SPA and 4PA, accuracy improves since 
deviation decreases to 5.1 10~^ at C3 = 1- 



A. Samples with strong linear absorption 

Literature values for the linear-optical absorption coef- 
ficient ai in the near-infrared region vary between 0.1 and 
100 cm~^ for PTS,[25] the example given above, giving 
characteristic lengths as short as Li = 100 jam. Linear 
absorption dominates over NLA processes at sufficiently 
low values of /q and, therefore, such an effect cannot 
be neglected in this regime. In this sense, a solution of 
Eq. ^ may be provided for / = 1, that is, when linear 
absorption and only a single nonlinear process is active, 
which is written as UM, uM 



Here, 



Ceffr7 



exp(-7^^lCm) 

+ - l)Ceffm 

1 - exp [-(m - l)r]miCm] 
(m - l)r]mi 



(IS) 



(14) 



stands for the effective axial coordinate. Importantly, 
the product r]miCm (— Ci) is independent of Lm- Also, 
Ceffm and Eq. ([13]) lead to Cm and Eq. (|4]), respectively, 
if Vmi = 0. 

Analyticity of the solution ([13]) suggests that an ap- 
proach of / may be given for the complete equation 



dcJ=-VmJ-i'^-Vmll\ 



(15) 
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FIG. 3: Dependence of / at different depths when mPA, /PA, 
and Unear absorption are found simultaneously at r/mi = 1 
and r]mi = 1- Again, the numerical estimation of / is plotted 
in solid lines, and the approximation given in ()17|) is drawn 
in dotted Unes. 



beam is required. Also, Gaussian beams are found 
more attractive in theoretical studies because the three- 
dimensional distribution of the field is an analytic so- 
lution of the paraxial wave equation. [28] In this study 
we assume a sufficiently narrow bandwidth of the pulsed 
wave in order to neglect spatiotemporal couplings associ- 
ated with ultrafast beam propag ation.[29',^] For conve- 
nience, we use normalized spatial coordinates: the radial 
coordinate r involves a normalization with respect to the 
beam waist wq and the axial coordinate z is given in 
units of the Rayleigh range zr = kwQ/2 being k = uo/c 
the wavenumber in air. If the waist plane of the Gaussian 
beam is found at z = (neglecting focal shifts, [31] this is 
the focal plane of the lens), the instantaneous intensity 
at the input plane of the scanned sample yields 



/o(r, z, t) = E(t)w~'^ exp (-2rV^^ 



(18) 



being 



which accounts for linear absorption and two different 
NLA processes, simultaneously. Accordingly, we may 
follow the procedure given above identifying the terms 
in the parabolic series expansion of / about the ori- 
gin, = 0, derived here from Eq. ([13]) and also from 
Eq. (p!5|) . Using over Eq. (p!3|) the normalization of the 
axial coordinate Cm L/Lmi in terms of the average 
length of Eq. ([8|), and also substituting the number m 
of photons by the average ffi of (pT]) , lead to severe sim- 
plifications. The required approach is straightforwardly 
attained if we additionally substitute rjmi = Lm/Li by 
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giving 



exp {-rjmlCm) 



-^l + (m-l)(l + r/™OC 



(16) 
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The effective axial coordinate Ceffm is also conveniently 
transformed from (p!4|) by using the average ffi instead of 
the integer m. 

Validity of Eq. (pT|) is examined in Fig. [3] for balanced 
absorption processes, i.e., rjrni = 1 and rjrni = 1- Nu- 
merical simulations are also depicted, and an excellent 
agreement with Eq. ([T7j) is remarkable. We point out 
that the approach (pT|) provides the exact value of / in 
the limits r]mi = 0, r]mi = 0, rjmi oo, and rjmi oo. 



III. Z-SCAN SIGNATURES WITH GAUSSIAN 
BEAMS 

From an experimental point of view, nonlinear sam- 
ples are commonly proved with paraxial Gaussian 
beams focused with low- numerical- aperture lenses. Top- 
hat beams and trimmed Airy beams are claimed to 
present some advantages over Gaussian beams for Z- 
scan experiments [27]; however, beam shaping of the laser 



w{z) = vTT" 



(19) 



the Gaussian width, and E{t) the intensity of the pulse 
envelope at focus, {r^z) = (0,0). The temporal coordi- 
nate t supports a normalization over the pulse duration 
such that 



dtE{t) = 1, 



(20) 



where Eq is the maximum instantaneous intensity. For 
the sake of convenience, from hereon we consider pulse 
waveforms with a Gaussian profile. In this case we write 
E{t) = Eoexp{-7rt^) so that its FWHM yields Ar = 
2Y^ln2/7r 0.94. We point out that our analysis may be 
extended to other waveforms straightforwardly. 

In the open-aperture mode of Z-scan measurements, 
light transmitted through the sample is fully collected. 
Mismatching of the linear index of refraction of the sam- 
ple and air induces the incident beam is partially re- 
flected, an effect that may be treated on the experimental 
data and is therefore ignored here. The output intensity 
/ = IqI is obtained from Eqs. (pT|) and (p!8|) , which varies 
in space and time. For a given pulse, the coordinate Cm 
reaches a maximum value 



Am = OLmLE'^ ^ 



(21) 



at focus. Also the coupling parameter satisfies r]mi < 
Ai/Ajji (if m < I). Upon recording of the integrated 
intensity at the exit plane of the scanned sample. 



/oo poo 
dt I 2iTrdrI(r,z,t), 
-oo 



(22) 



changes are found in the vininity of the Gaussian focus, 
|z| < 1. In the far-field, \z\ ^ 1, the output intensity 
/ ^ /o exp (-Ai) so that Eq. ([22]) yields 



Poo = lim P{z) = -Eo exp (-Ai) . 

z — >oo Z 



(23) 
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FIG. 4: Z-scan signatures for 2PA and SPA admixes (left) and simultaneous SPA and 4PA processes (right). Dominant mPA 
with Am = 1 (and A; = 0) is drawn in solid line. Other values of the ordered pair (Am, A^) are represented in different curve 
styles. 



Commonly, Z-scan signatures refer to traces of the ra- 
tio T = P/Poo, which is interpreted as the normalized 
transmittance of the sample as being scanned axially. 

In order to simplify our discussion remaining in 
Sec, mil let us neglect linear absorption by using the limit 
Ai 0. In Fig. m we depict Z-scan signatures when the 
sample experiences 2PA and SPA simultaneously (subfig- 
ure on the left), and when SPA and 4PA are present (on 
the right). The strength of each NLA process is char- 
acterized by the corresponding parameters Am and A^ 
Numerical simulations of the transmittance T are per- 
formed using the approximated equation (p!2]) : we have 
observed by direct comparison with numerical computa- 
tion of Eq. ([1]) that these traces are of extremely-high 
accuracy for low and intermediate values of A^ and A^ 
Z-scan responses are symmetric even functions with re- 
spect to the origin, z = 0. As expected, NLA is maximum 
at the Gaussian focus. The V-shaped response T shows 
a distinct tail decay and valley value, a fact that may be 
employed to determine NLA coefficients of the sample 
from experimental data. However, a procedure may be 
given to overcome the inherent lack of precision in the 
data fit when simultaneous NLA are present, specially as 
the number of absorbed photons in each nonlinear pro- 
cess increases. 

We may take advantage of Eqs. ([7j) and ([9]) in order to 
obtain analytical approaches of Z-scan signatures. Af- 
ter denormalizing /^^^ conveniently and inserting into 
Eq. ([22|) , we finally have a first-order T response. 



= 1 



,.2(m- 



,2(1-1)^ 



(24) 



Using /^^^ instead, a second-order approach of the Z-scan 
signature yields 



2^(2) ^.(i) 



+ 



(/ + m)A^A, 



2^2(m+Z-2) ^(r^i + /-l)3 



(25) 



2^4(/-i) ^(2/ - 1)3 2^4(m-i) ^(2m- 1)3 



Again, in Eq. (p!8|) we have considered a bell-shaped 
Gaussian pulse. Following a similar approach, Eric W. 
Van Stryland et al.[S] have provided an infinite series 
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FIG. 5: T(z) response evaluated numerically (solid curves) 
and from the linear approach (|24|) . T^'^\ in dashed lines and 
the quadratic approximation (|25)) . T^'^\ in dotted lines. 



for a solitary 2PA process (m = 2 and Ai = 0). For 
simultaneous mPA and /PA we might give also an ex- 
pansion of infinite number of constituents; however, the 
most relevant terms are provided in T^^^. We point out 
that top-hat pulses (and of some other shapes) yield also 
analytical expressions for the coefficients of the series ex- 
pansion ([25|) . 

High accuracy of T*^^^ is observed at A^, A^ < 10~^; 
this fact is expected also from numerical computations 
depicted in Fig. [H In these cases, however, the collected 
signal slightly deviates from the unit transmittance asso- 
ciated with the response in the linear regime. This effect 
is noticeable at increasing integers m and / where NLA 
is highly localized in space (around the focus) and time 
(about the local time t = 0). 



Values of A^ and Ai commonly exceed 10 ~ in empir- 
ical demonstrations of Z-scans, so we may test the per- 
formance of Eqs. ([24|) and ([25|) in those practical cases. 
We tentatively assume that deviations of the normalized 
transmittance T induced by self-focusing are negligible. 
In Fig. [5] we plot Z-scan signatures for mixed mPA and 
IFA of different orders, however keeping a balanced peak 
depth A^ = A^. A modest approach of T*^^) and T*^^^ 
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FIG. 6: Asymptotic behaviour of the absorptance A versus 
the Gaussian width w, drawn in sohd hues. Balanced depths 
Am = = 1/2 have been used in the simulations. Dashed 
lines correspond to asymptotes (|27)) . 



is observed for a value = 2 10~^; discrepancies are 
significant only at the Gaussian focus and nearest axial 
points. When A^ reaches the unity, the series expan- 
sion T^*^^ manifests a strongly-oscillatory convergence at 
increasing specifically in the focal region of the Gaus- 
sian beam. In all cases, however, an excellent agreement 
of T(i) and T(2) responses is found along the signature 
tails. 



IV. DETERMINATION OF THE NLA 
COEFFICIENTS 

In spectroscopic measurements, absorption coefficients 
are determined after data processing by fitting the re- 
sponse T within the region of interest; commonly this 
is the focal region of the Gaussian beam where the sig- 
nal changes stronger. The detected transmittance valley 
has a shape and a value at the bottom that may help to 
determine NLA coefficients; however, this method has in- 
trinsically poor sensibility for two or more unknowns and 
may lead to ill estimates. On other hand, far-field decay 
of Z-scan signatures, which is delineated accurately in 
Eqs. and ([25|) . reveal distinctively the NLA orders 
of the absorptions involved and the values of their respec- 
tive coefficients. As a consequence, data fitting may be 
carried out alternatively onto points of the trace out of 
focus, assuming that weak disturbances of the transmit- 
tance may be discerned in spite of fluctuations. 

Instead of using the normalized transmittance T expe- 
rienced in the sample as being scanned along the focal 
region, for convenience, our theoretical analysis is sim- 
plified in terms of the absorptance 



Aiz) 



T{z) 



(26) 



produced in the nonlinear material. Specifically, the 
curve A{z) traces a centered peak of symmetric vanishing 
tails. In accordance with ([24|) and ([25]) , and assuming 
that m < the absorptance tails mainly fall in direct 
proportion to w raised to the power of — 2(m — 1). For 
illustration, in Fig. [6] we show a log- log plot of the peak 



FIG. 7: Absorptance A (soUd lines) and its asymptotic curve 
(dashed lines) for opaque samples (Ai ^ 0) in the linear 
regime. Values Am = A; = 1/2 are employed in the sim- 
ulations for: (a) m = 2 and / — 3, and (b) m = 3 and / — 4. 



A versus the Gaussian width w. In the far field, w ^ 1^ 
the absorptance draws a straight line 



log A = — 2(m — 1) log w + log 



(27) 



given in the slope-intercept form. We point out that dis- 
crepancies of A and its asymptote ([27|) are not significant 
even at il' > 1 for small and moderate values of A^ in com- 
parison with Am- This asymptote has a negative slope 
of absolute value 2(m — 1), thus the integer m may be 
revealed unambiguously by linear curve fitting to exper- 
imental data. The intercept of the function ([27|) found 

w = 1 displays m~^/^A^ and subsequently the peak 
depth may be estimated. Ultimately the mPA coefficient 
am is determined from Eq. ([2T]) if the sample width L 
and peak intensity Eq are known. 

Linear absorption is being omitted in the geometric 
analysis of the absorptance. In Fig. [71 we investigate the 
asymptotic behaviour of A beyond the Gaussian focus 
for different values of Ai = aiL. In comparison with the 
asymptote of Eq. ([27|) valid for Ai = 0, opacity of the 
sample (in the linear regime) leads to a parallel asymp- 
tote; thus its slope is still characteristic of the order m 
of the nonlinear ity. However, its intercept at c<; = 1 that 
allows an estimate of A^ is undervalued. This is not 
surprising since the parameter A^ stands for the maxi- 
mum value of Cm, which is a coordinate that leaves out 
of consideration linear absorption. Otherwise, we have 
suggested above the use of the effective axial coordinate 
Ceffm given in Eq. ([M]) instead. In accordance we define 
the effective parameter Aeffm = Qm^m as the maximum 
value of Ceffm, where 



1 -exp [-(m- l)Ai] 
(m- l)Ai 



(28) 



denotes the rate decrease of the estimate A^, and there- 
fore the estimate decrease of a^n, in the presence of linear 
absorption. Fig. [7] also shows the corresponding asymp- 
totes of Eq. ([27|) when A^ is substituted by Aeffm, which 
is in agreement with our analysis. 

The steps developed above may be applied trivially for 
samples experiencing only one NLA process (A^ = 0). 
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FIG. 8: Signal A- Am/w^^"^-^^ 
tote (|3Qp in dashed lines for Ar 



,3/2 
A, 



Signal A — Am/w'^^^ ^^m'^^^ in solid lines and asymp- 

1/2. 



To the knowledge of the author, this is the first time a 
procedure is proposed to derive am asymptotically from 
a single Z-scan trace. Interestingly, in Refs. @, [Til, and 
[32I , alternate asymptotic analyses are proposed for single 
NLAs based on absorptances measured at the Gaussian 
focus for different peak intensities Eq of the input beam. 



Using w = 1 and amLE^ , Eq. ([27j) leads to 



log A = (m — 1) log£^o + log 



am.L 



(29) 



Once again, the order m and the coefficient am may be 
determined geometrically in a log-log representation of 
A versus £^0 from the slope and intercept of the asymp- 
tote ([29|) . respectively. If linear absorption cannot be 
neglected, we would find a multiply factor Qm to correct 
the value of am- In the case that two NLA processes are 
manifested simultaneously, the conclusions given above 
are already valid assuming that m is the lowest order of 
the nonlinearity. 

Following with light-matter interactions involving an 
admix of two different NLA processes, we emphasize that 
accurate measurements leading to a precise Z-scan trace 
should contain sufficient information to extract not only 
the order m and the NLA coefficient am but, additionally, 
estimates of the parameters / and a^ For that purpose, 
we propose to substract the foremost term of T^^^ (apart 
from unity), A^/i(;^*^^~^^A/m^, yet considered, from the 
absorptance peak A. The tail of the resultant signal has a 
characteristic decay, faster than the removed part, which 
depends on the remaining NLA competing in the sample. 
In Fig. Owe illustrate this procedure with some numerical 
examples. In the far field, w ^ 00^ the signal exhibits an 
asymptote 



log 




-2{l - l)\ogw 



mA'^S2mhi 



(30) 



where S denotes the Kronecker's delta. In a log- log plot, 
the slope of (|3Q|) displays the order / of the second non- 
linear process and the intercept includes the value of 
the peak depth A^ In the particular case m = 2 and 




0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 

A A 



FIG. 9: Normalized peak absorptance Amax for different pa- 
rameters Am and A; and different NLA orders. 



/ = 3, however, the intercept incloses a term depending 
on the squared of A^. This is not a major problem since 
A2 is estimated beforehand. We point out that if also 
A2 = A3 = 1, which by the way are values in the bound- 
aries of validity of our approach, then both terms are 
equivalent in magnitude but having opposite signs and 
they cancel each other out. As a consequence, higher- 
order terms of the series expansion of T would dominate 
the asymptotic behavior of the signal. 

From the numerical simulations shown in Fig. [8] we 
are concerned about the extreme accuracy of the trans- 
mittance data, at least in the order of 10~^-10~^, re- 
quired to draw the trace of the signal in the asymptotic 
regime. Unfortunately, experimental measurements com- 
monly lacks of sufficient precision to apply this method. 
In these cases, data fit is suitably performed along seg- 
ments of the optical axis where absorptance is higher, for 
instance, at the Gaussian focus. At this point, addition- 
ally, the contribution of the higher-order NLA to the ab- 
sorptance peak, denoted as Amax^ is enhanced. Asymp- 
totic analyses of Amax — Am/Vrn^ in terms of ^o, in the 
way carried out in Eqs. ([29|) and ([3Q|) at w = 1^ are avail- 
able in order to estimate / and A^, but would demand 
also a strong agreement with the experimental data. 

Assuming that the order / of the nonlinear process is 
known, the peak depth Ai and thus the NLA coefficient 
ai is preferably determined by giving the best fit of Amax 
to data. In Fig. [9] we depict contour plots of the peak 
absorptance Amax for degenerate NLA of different orders 
in terms of the peak depths. The plots show a monotonic 
smooth variation allowing the estimation of A/, which 
represents the only unknown remaining in our problem, 
by means of a data fit. 

Theoretical evaluation of Amax displayed in Fig. [9] is 
carried out numerically since we cannot provide a sim- 
ple analytical expression of the absorptance dit w = 1. 
However a series expansion may be obtained from T^^^ 
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as 5* ^ oo, which yields 



{p + qV- ™ 

p,q=l ^ ^ 

Inserting Eq. ([5]) into Eq. ([22|) leads to 



(31) 



CM = 



(-1) 



9-1 



[qn - (q - l)]3/2 



T|[,^ _(,_!)], (32) 



s=0 



which are characteristic coefficients of individual NLA 
processes of order n = m and being g = 1, 2, ... a posi- 
tive integer. Additionally, an infinite number of coupled 
coefficients arise in (|3T]) , some of those of the lowest order 
being 



Cii(m, /) = 
C2i(m,/) = 
C3i(m, /) 
C22{m,l) 



m + / 



(m + /- 1)3/2 
Sm^ + 2m/ + /2 - 2m - / 



(33a) 
(33b) 



(2m + /- 2)3/2 

(4m^ + m/ + - 3m - (3m + / - 
(3m + / -3)3/2 

(Tm^ + 4m/ + 7/2 - 6m - 6/) (m + /^-^^)^ 



(2m + 2/ -3)3/2 



The property of symmetry Cpq{m^ I) = Cqpil^ m) allows a 
faster computation of different coefficients. As inferred in 
a previous discussion, as long as the peak depths ap- 
proach the unity, accurate estimation of Amax is expected 



to involve running terms of increasing index S = q'm 
the summation. 



V. CONCLUSIONS 

We have performed a geometrical analysis of the open 
aperture Z-scan trace when multiphoton absorption of 
different nature arise, simultaneously, in the material un- 
der inspection. Under appropriate conditions, a simple 
analytical expression may be found allowing an accurate 
description of the beam shaping induced by absorption. 
Parametrization of the nonlinear beam propagation is 
formulated in terms of a series expansion providing an 
analytical representation of the absorptance produced in 
the nonlinear medium as experimentally recorded in the 
traces. Unlike fitting Z-scan signatures directly, we pro- 
pose the evaluation of asymptotes arisen when absorp- 
tance are represented versus the width of the Gaussian 
beam at different scanning tracks. A fundamental draw- 
back of the procedure is the extreme sensitivity of data 
measurement required to draw the estimates out for the 
coefficient of the highest order. To overcome this ma- 
jor inconvenient, the peak value of the absorptance is 
suggested to play on the determination of the remaining 
NLA coefficient. 
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